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ABSTRACT 

We use high-resolution (0.5 arcsec) CO(2-l) observations performed with the Atacama Large Millimetre/submillimetre Array to trace 
the kinematics of the molecular gas in the Seyfert 2 galaxy IC 5063. The data reveal that the kinematics of the gas is very complex. A 
fast outflow of molecular gas extends along the entire radio jet (~ 1 kpc), with the highest outflow velocities about 0.5 kpc from the 
nucleus, at the location of the brighter hot-spot in the W lobe. The ALMA data show that a massive, fast outflow with velocities up to 
650 km s“’ of cold molecular gas is present, in addition to one detected earlier in warm H 2 , H i and ionised gas. All phases of the gas 
outflow show similar kinematics. IC 5063 appears to be one of the best examples of the multi-phase nature of AGN-driven outflows. 
Both the central AGN and the radio jet could energetically drive the outflow. However, the characteristics of the outflowing gas point 
to the radio jet being the main driver. This is an important result, because IC 5063, although one of the most powerful Seyfert galaxies, 
is a relatively weak radio source (P, 4 qhz = 3 x 10^^ W Hz“*). 

All the observed characteristics can be described by a scenario of a radio plasma jet expanding into a clumpy medium, interacting 
directly with the clouds and inflating a cocoon that drives a lateral outflow into the interstellar medium. This model is consistent with 
results obtained by recent simulations such as those of Wagner et ah. A stronger, direct interaction between the jet and a gas cloud is 
present at the location of the brighter W lobe. This interaction may also be responsible for the asymmetry in the radio brightness of 
the two lobes. 

Even assuming the most conservative values for the conversion factor CO-to-H 2 , we find that the mass of the outflowing gas is between 
1.9 and 4.8 x 10^ Mq, of which between 0.5 and 1.3 x 10^ Mq is associated with the fast outflow at the location of the W lobe. These 
amounts are much larger than those of the outflow of warm gas (molecular and ionized) and somewhat larger than of the H i outflow. 
This suggests that most of the observed cold molecular outflow is due to fast cooling after being shocked. This gas is the end product 
of the cooling process, although some of it could be the result of only partly dissociated clouds. Our CO observations demonstrate 
that fast outflows of substantial masses of molecular gas can be driven by relativistic jets, although in the case of IC 5063 the outflows 
are not fast enough to remove significant amounts of gas from the galaxy and the effects are limited to the central ~ 0.5 kpc from the 
centre. 
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1. Introduction 

The discovery of the large variety of structures and phenomena 
traced by cold gas in early-type galaxies has made these objects 
much more complex than historically perceived. Although not a 
compl ete su rprise (see e.g. t he wo rk of ivan Driel & van WoerdenI 
i 199 ill and IWiklind et al.l (Il99-5h . the plethora of new results, 
made possible by the technical improvements in the available 
and new observing facilities, has opened new and important 
avenues for the study of formation and evolution of early- 
type galaxies us ing this component of their inter-stellar medium 


l2013UAlataloetaLll2013: 

Serra et 31.12014. 

le formation his- 
nexpected char- 
for tracing ex- 
ust AGN-driven 
glio et alJ 2 OIOI: 

In addition to being an effective tracer of tl 
tory of galaxies, the cold gas has revealed u 
acteristics that have allowed to use it also 
treme energetic phenomena like massive and 
gas outflows (Morganti et al. 2005 alibi: Peru 

iDasvra & Combesll2012HMorganti et alil2013al 

btlMahonv et alJ 


I 2 OI 3 I : lMorgantill20l4 and refs therein). This finding has chal¬ 
lenged our ideas of how the energy re leased by an AGN may 
interact with its surroundin gs (see, e.g.. lTadhunteill2()(MlFabianl 
I 2 OI 3 I : l^mbes et ^l2014l for reviews on the topic). Cold gas 
has been found to represent the do minant mass compone nt 
of the outflows observed (s e e, e.g., Morganti et alJ I2005alf9: 
Kaneky &Chengaluil 200^ Holke^^ 20091 iFeruglioet^p 
20101: lAlataloet al.M201iriDasvra&Combesl I 2 OI 1 I 120121: 


Morganti et alJl201^bl : ICicone et alJ2014l:lGarcfa-Burillo et al.l 

2014 for some examples). This makes them particularly inter¬ 
esting for their potential impact on the ISM of the host galaxy. In 
the effort of understanding the mechanisms that can regulate the 
black-hole growth and the quenching of star-formation in early- 
type galaxies, the discovery, in a growing number of objects, of 
these AGN-driven massive outflows represents an important step 
toward building a more complete and realistic view of how feed¬ 
back may work. Thanks to detailed studies of these phenomena 
in nearby objects, better constraints (i.e. distribution and kine¬ 
matics of the outflowing gas, energy injected into the ISM, etc.) 
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Fig. 1. Total intensity image representing the distribution of the CO(2-l) in IC 5063 and showing the striking brightness contrast between the 
inner, bright CO and the fainter outer disk. The details of the bright inner region can be seen in Fig.|2] The resolution of this image is 0.54 x 0.45 
arcsec (PA = -31.7°). This image is not corrected for the primary beam of ALMA. 


can be provided to theoretical models of galaxy formation that 
now routinely include feedback effects. 

Different mechanisms have been proposed to accelerate the 
gas. Although (broad) winds from the accretion disk, launched 
by radiation pressure, interacting and shocking the surrounding 
medium is often f avoured (e.g. IZubovas & King! 1201 2l 120141: 
ICosta et ^ l2014t ICicone et al.l l2014h . the role of the radio 
plasma has also gained interest. Indeed, this mechanism can 
be particularly relevant in early-type galax ies where up to 3 0% 
of the high-mass g alaxies are radio loud dBe st et al. 2005). A 
number of studies ( Birzan et"^ 120081: ICavagnolo et al.l 20101: 
iMcNamara & Nulsenll2012|) have shown an high efficiency of 
the coupling between the radio plasma and the surrounding 
ISM/IGM, and mechanical power of the radio jets exceeding the 
synchrotron power. These results suggest that this mechanism 
can be relevant also for relati vely weak radio sources (see, e.g., 
NGC 1266 and NGC 1433. lAlatalo'eFaDl201 ll: ICombes et al'l 
120131 respectively). 

Expanding the number of known outflows of atomic (H i) 
and molecular gas in order to provide bett er statistics is the f o¬ 
cus of a number of r ecent studies (see, e.g.. lGereb et al.l2014allbt 
ICicone et al.l l2014l) . However, the other important direction 
where progress needs to be made is to spatially resolve the out¬ 
flows and image in detail the distribution, kinematics and the 
physical properties of the gas. In the case of 4C 12.50, VLBI ob¬ 
servations have, for the first time, allowed this t o be established 
for the Hi component (iMorganti et akl I^13ah . For the warm 
molecular gas, the most spectacular and convinc ing example of a 
j et-dri ven outflow has been found in IC 5063 bv lTadhunter et al.l 
(l2014l) . In this object, the warm component of the molecular gas 
has been identified at the location of the interaction between the 
radio jet and the ISM, 0.5 kpc from the central A GN, where an 
Hi outflow was also found (lOosterloo et al.ll2000l) . AGN-driven 
outflows of c old molecular gas have been so far only barely 
resolyed (e.g. iFeruglio et al.l l20f^ ICicone et al.l 1201^ not al¬ 
ways allowing a detailed modelling of the outfl owing gas and its 
conditions. The possible exceptions are M51 dMatsushita et al.1 


200711201 5h and NGC 4258 (iMartin et al.l 1 19891: iKrause et al.l 
20071) . The situation is now chang ing with the adyent of ALM A. 
The detailed case of NGC 1068 dGarcfa-Burillo et ^l2014l) is 
an example of the dramatic improyement that ALMA allows to 
make. 

All together, these studies are the first to proyide details on 
the physical parameters and the conditions of the gas in the out¬ 
flows. One of these objects, the Seyfert 2 galaxy IC 5063, is the 
target of the ALMA obseryations presented here showing, for 
the first time, the complexity of the kinematics of the cold gas 
under the effect of the radio jets. 


2. A case study: IC 5063 


IC 5063 is a Seyfert 2 object hosted by an early-type galaxy 
with a prominent dust lane which is seen edge on. Together with 
NGC 1068, it is one of the most radio bright Seyfert objects 
with a poweiQ of Pi. 4 ghz = 3 x 10^^ W Hz ' dTadhunter et al.l 
l2014l) . In the radio continuum, IC 5063 shows a triple struc¬ 
ture of about 4 arcsec in size (about 0.93 kpc; see lMorganti et al.1 
d 19981) ) aligned with the dust-lane. Hi emission (8.4 x 10® Mq) 
is traced in a regularly rotating warped disk with radius ~ 2 ar- 
cmin (about 28 kpc) and which is associated with the system 
of dust lanes. The likely counterpart of the Hi disk is also de¬ 
tected in CO as orig inally obseryed with SEST in CO(l-O) by 
[WMindet^ (1^^) and more recently with APEX in CO(2-l) 
by iMorganti et al.l (I 20 1 3bl) . 

IC 5063 was the first object where a fast, AGN d riyen mas- 
siye outflow of H i was detected (IMorganti et al.ll99^ . This out¬ 
flow was found to be located agai nst the brightest radio hot-spot, 
about 0.5 kpc from the radio core (lOosterloo et al.l2000t) and this 
has been interpreted as eyidence that the outflow is driyen by the 
interaction of the radio jet with the surrounding ISM where the 


' Assuming a Hubble constant H„= 70 km s“* Mpc^' and Oa = 0.7 
and Am = 0.3, we use an angular scale distance to IC 5063 of 47.9 Mpc, 
implying a scale of 1 arcsec = 232 pc. 
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Fig. 2. Contours of the 230 GHz continuum emission superposed to the 
central region of the total intensity of the CO(2-l) illustrating the spatial 
correlation between the two. Contour levels are 90, 180, 360, 720, 1440, 
2880, 5760 and 15120 ^liy beam-*. 


jet is moving in the plane of the large-scale gas disk of IC 5063. 
This hypothesis has been further strengthened by observations 
(using ISAAC on the VLT) of warm molecular gas (H 2 at 2.2 
yum) which have, for the first time, shown the presence of out- 
flowing molecular gas co-spatial with the bright radio hot-spot 
dTadhunter et al.ll2Q14l) . A signature of outflow of cold molecu- 
lar gas has been found using CO(2-l) observations with APEX 
(iMorganti et al.ll^l3bl) . The limited spatial resolution of these 
observations did not allow to determine where the outflow is oc¬ 
curring in relation to the radio structure and the main goal of 
the ALMA data presented in this paper is to remedy this situa¬ 
tion. Th e data presented here complement the in formation about 
ionize d (iMorganti et al.l2007h. warm molecular ( Tadhunteret_^ 
2014 ) and cold atomic (iMorgantiet 119981: lOosterloo et alJ 
200 c ). making IC 5063 the best case of a multi-phase AGN- 
driven outflow studied so far. 

The sensitivity and spatial resolution of the new ALMA ob¬ 
servations allow us to show for the first time the detailed distri¬ 
bution and kinematics of the cold molecular fast outflow, provid¬ 
ing a major step forward in the understanding of the interaction 
between the energy released by the AGN and the surrounding 
medium. 


Fig. 3. Integrated position-velocity map taken along the major axis of 
IC 5063 where the data were integrated over a length of 2 arcsec in 
the direction perpendicular to the major axis. T he white line gives the 
rotation curve we derived from the photometry of lKulkarni et all (1199 8h 
and indicates the expected kinematics of gas following regular rotation. 
Contour levels are -3, -1.5, 1.5, 3, 6, 12 and 24 a. Positive radii are east 
of the centre. 


data were subsequently exported to MIRIAE0 (ISault et al.lll995h 
which was used to perform additional self-calibration which im¬ 
proved the quality of the images significantly. All further reduc¬ 
tion steps (continuum subtraction, mapping/cleaning) were also 
done in MIRIAD. 

The data cubes were made using various Briggs weightings 
in order to explore the optimum for imaging. However, the im¬ 
ages did not show any major differences (owing to the good uv 
coverage of ALMA). Here we present the final data cubes ob¬ 
tained using Briggs robust = 0. The line cube was produced us¬ 
ing a velocity resolution, after Hanning smoothing, of 20 km s *. 
At this velocity resolution, the r.m.s. noise per channel is ~ 0.3 
mJy beam '. The large array configuration used is providing a 
resolution of the line data of 0.54 x 0.45 arcsec (PA = -31.7°), 
nicely allowing the line emission to be resolved in relation to the 
radio continuum emission. 
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3. The ALMA data 

IC 5063 was observed during Cycle 1 on the 27 May 2014 with 
the Atacama Large Millimeter/submillimeter Array (ALMA) 
using Band 6, covering simultaneously the CO(2-l) centered 
on 227.968 GHz and two additional base-bands (at 232 and 
247 GHz) for imaging of the continuum. The observations were 
pointed at the nucleus of IC 5063, with a single-pointing held 
of view (LoV) of ~ 27" and making use of the correlator in 
Lrequency Division Mode with a 1.875 GHz bandwidth (cor¬ 
responding to 2466 km s ' at the frequency of CO(2-l)). The 
frequency resolution was 976.5 kHz (1.28 kms '), but channels 
were combined when making the image cube to provide a bet¬ 
ter match to the observed line widths. The observation made use 
of 31 antennas with a maximum baseline of 1.5 km. The dura¬ 
tion of the observations was 1 h which was required to reach the 
sensitivity to locate and trace the outflow of m olecular gas. The 
initial calibration was done in CASA0 (v4.2 iMcMulline et al.l 
l2007h using the ALMA reduction scripts. These calibrated uv 


In Lig. [T]the total intensity of the CO(2-l) is shown. The 
single pointing ALMA observations detects molecular gas up to 
a radius of ~ 10", covering therefore only the very inner part 
(r < 2 kpc) of the H i disk observed with ATCA (IMorganti et al.l 
Il998h . The bulk of the emission covers the same range of veloc¬ 
ities (betw een 3140 to 3630 kms~ ') as traced by the APEX ob¬ 
servations (IMorganti et al.ll^l3bll . However, the wider band of 
ALMA allows to trace a much broader component of blueshifted 
gas, going down to velocities of ~ 2750 kms^*, i.e. up to ~ 650 
km s“' blueshifted compared to the systemic velocity. 

A continuum image was also produced from the 2 additional 
basebands using uniform weighting. The r.m.s. noise of the con¬ 
tinuum image is 30 luJy bearn^*. The restoring beam 0.45 x 0.42 
arcsec (PA= -30.4°). Ligure |2] shows the comparison between 
the mm continuum and the molecular gas in the inner region. 

Linally, the CS(5-4) line at 242.205 GHz was included in a 
fourth base-band, but no emission from this spectral line was 
detected. The noise level in this cube, not quite uniform across 
the band, ranges between 0.3 and 0.4 mJy beam”'. 


^ http://casa.nrao.edu/ 


^ http://www.atnf.csiro.au/computing/software/miriad/ 
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4. Results 

The main result of our observations is the detailed view of the 
complexity of the distribution and kinematics of the gas in the 
central regions. The data clearly demonstrate that, apart from 
the regularly rotating large-scale disk, the distribution and kine¬ 
matics of the gas in the inner regions is strongly affected by the 
radio jet. 

Figure [1] shows that the distribution of the CO(2-l) follows 
two structures: a high-surface-brightness inner structure of about 
the size of the radio continuum and a lower-surface-brightness, 
larger disk. This large disk is likely the molecular counterpart of 
the Hi disk known to be present in this galaxy (see Sec. 14.3b . 
Some of the gas in this disk appears to be distributed along arm¬ 
like features. However, the relatively edge-on orientation of IC 
5063 makes it difficult to identify and characterise these struc¬ 
tures well. 

The bright and elongated structure observed in the central 
regions shows a striking correspondence with the radio contin¬ 
uum, providing the first clear signature of the impact of the radio 
plasma on the ISM. Figure|2]illustrates how the CO bright emis¬ 
sion wraps around the radio continuum, with the brightest CO 
component coincident with the brighter radio lobe on the west¬ 
ern side. 

Below we describe kinematics and other properties of these 
CO components. 




2400 2900 3400 3900 



over the western radio lobe of IC 5053 are compared with th e spatially 
integrated Hi 21-cm absorption ifrom lMorganti et'^ll998h . The plot 
shows that all the components have gas reaching velocities of -650 
km s“' or beyond. 


4.1. A complex, extended outflow of molecular gas 

The kinematics of the gas provides a striking view of what is 
happening in the central region. In Fig. |3] we show a position- 
velocity ipv) diagram of the CO gas taken along the major axis 
of the galaxy (which coincides, to within a few degrees, with the 
orientation of the radio continuum source). Given the relatively 
edge-on orientation of the system, a position-velocity diagram 
gives a better view of the kinematics of the gas compared to a 
velocity field. To make this pu diagram, we have integrated the 
data over a length of 2" in the direction perpendicular to the 
major axis. Therefore this map shows all the gas in the vicin¬ 
ity of the radio source, but some of the gas associated with the 
large-scale disk is not included. To indicate what can be expected 
for the kinematics of regularly rotating gas, we have overplotted 
a rotation curve base d on the HST photometry of IC 5063 of 
iKulkarni et al.l (Il998h (see below). 

Figure[3]clearly illustrates the strikingly complex kinematics 
of the gas in the region co-spatial with the radio plasma. Beyond 
a radius of 2" (i.e. ~ 460 pc), the gas follows the expected rota¬ 
tion pattern fairly well and this gas corresponds to the large-scale 
disk, unaffected by the radio plasma. If gas with such complex 
kinematics would have been present outside the region of the ra¬ 
dio emission, the sensitivity of our observations would be able to 
see. Instead, the change between complex to regular kinematics 
is quite abrupt. However, inside this radius, most of the gas has 
velocities quite different from what is expected based on the ro¬ 
tation curve. In fact, a significant fraction of the gas is found in 
the so-called "forbidden quadrants" of the pu diagram, i.e. those 
parts where co-rotating gas cannot appear. Examples of this are 

the gas at r-2" and with velocities around 3500 km s“* and, 

similarly, the gas at r ~ 1" and V ~ 3200-3300 kms^'. In ad¬ 
dition, for -2" < r < 0" there is a substantial amount gas with 
velocities that are well outside the range of expected velocities, 
indicating that the kinematics of this gas is not determined by 
gravity alone. This gas is most likely related to t he fast, outflow- 
ing Hi gas that was discovered in absorption bv iMorganti et al.l 


(Il998l) . As we will show below (Sec. I5.2l i. the distribution and 
kinematics of the molecular gas with anomalous velocities un¬ 
equivocally point to the presence of outflowing gas and can be 
understood in terms of a lateral outflow driven by the radio jet. 


In further ag reement with earlier observations, also of 
other gas phases (iMorganti et al.l 20071: Tadhunter et al.l 120141: 


iMorganti et aklfToOsT Oosterloo et al.f 2000h . the most extreme 

velocities are confirmed to be present at the location of the bright 
W radio lobe, about 0.5 kpc from the nucleus. This, again, is 
a strong indication that the most likely origin of these extreme 
kinematics must be the jet. The most extreme velocities of the 
anomalous gas are blueshifted by more than 500 km s ' and 
redshifted by more than 400 km s ' w.r.t. the regularly rotat¬ 
ing gas at that location. The velocity profiles of the different 
gas components integrated over the region of the western hot¬ 
spot (except for the H i that shows the profile integrated over the 
full radio continuum region although the most blue shifted HI 
absorption is most likely to occur only against the continuum 
of the bright W hot spot) are compared in Fig. |4] In all cases 
(HI absorption, the ionized gas and the warm and cold molecu¬ 
lar gas), the bluesh ifted component of the gas reaches velocities 
of ~ -700 km s ' (IMorganti et al.l 1 199^ 120071: fTadhunter et al.l 
1201 4 respectively). 


The rotation velocity of IC 5063 is about 240 km s ' which 
means that the escape velocity in the inner regio ns is likely to 
be above 550 km s ' (iBinnev & Tremain3l2008l) . We observe 
basically no gas with such extreme velocities (see Fig. Ell. This 
means that, allowing for projection effects, at most a very small 
fraction of the anomalous gas will leave IC 5063 and that any 
possible effect on future star formation will have to be through an 
increased energy input in the ISM and not through gas removal. 
A si milar situation has b een observed to be the case in NGC 
1266 (lAlatalo et al.ll2015l) . 
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4.2. The bright inner CO(2-1) distribution 

A characteristic of the CO(2-l) distribution that clearly appears 
from Fig. [T] as well as from the position-velocity diagram in 
Fig. [3] is that, although the large-scale, regular rotating disk is 
clearly seen, the emission from the region around the radio jet is 
much brighter in CO(2-l). Two possibilities can be suggested to 
explain this. 


The bright inner CO (J=2-l) emission in IC 5063 can 
be excited by a variety of physical processes, such as 


shocks (e.g. ' 

Tollenbach & McKee 19891). cosmic rays 

(e.g 

Sternberg et al 

19981: Ferland et al.ll2009l IMeiierink et al.l 

2011 

Bayet et alJ 201 Ih. mechanical heating (iLoenen et alJ 

200 a 

Meiierink et alJ 20111) and photodissociation (including X- 

rays; e.g. Tielens & Hollenbachl 119851; iMaloney et al.l Il996t 

Meiierink et alJ 2007). Giyen that we haye measured only one 


molecular line (CO J=2-l) at sufficient spatial resolution to re¬ 
solve the molecular gas in this region we not able to discriminate 
between these different processes. 


The difference between the inner and outer regions may 
reflect differences in the temperature and/or excitation of the 
molecular gas in the two regions due to the influence of the 
jet/cloud interaction. As discussed by iBolatto et al.l (l2013h for 
starbursts and ULIRGs, to first order, higher gas temperatures 
yield brighter CO emission, decreasing Xco for the regions af¬ 
fected by the interaction. However, the increase i n density may 
partly compensate for this effect, as pointed out by lBolatto et al.l 
(1201!^ . Therefore a certain leyel of compensation may occur, 
lessening the impact of enyironment on the conyersion factor. 
Without a detailed study of the physical conditions of the molec¬ 
ular gas, e.g., using the full CO ladder, it is not possible to fully 
quantify these effects. 

On the other hand, the kinematics of the gas is clearly in¬ 
dicating that the effect of the radio jet is releyant for the distri¬ 
bution of the gas in the inner regions. The cocoon of shocked 
material created by the jet moying through the ISM may be able 
to push the gas and moye it outward. This explanation would 
be consistent with what i s proposed by theoretical models (e.g. 

I Wagner & Bicknelill201 ih and we will consider such a scenario 
in more detail in Sec. 15.21 

This scenario also implies that shocks will likely play a 
prominent role in the heating and excitation of the molecular 
gas. In addition the outward moying jet and associated shocks 
will boost the ionizing flux of cosmic rays (CR). Seyeral theo¬ 
retical studies haye shown that high CR rates can significantly 
modify the chemis t ry and conditions of the molecular gas (e.g. 

iFerland et alj|2009l : iBayet et ^120111 : iMeiierink et al.ll20Tlh . 

4.3. The reguiariy rotating CO disk 

Outside the region of the radio continuum emission, the CO(2- 
1) shows regular kinematics. The characteristics are of a rotat¬ 
ing disk with the western side approaching. This structure most 
likely represents the molecular counterpart of the regular disk 
already obseryed in Hi and in ionized gas (iColina et al.1ll99ll : 
iMorganti et aklflO^ l2007h . Indeed, the rotation yelocities at 
large radii of this molecular gas are yery similar to those of the 
H I, with an amplitude of the rotation of ~246 km s '. Howeyer, 
giyen the limited field-of-yiew of the primary beam of ALMA, 
we cannot determine the full extent of the molecular component 
of the large-scale disk. Furthermore, the comparison with the ye- 
locity profiles of the H i in the inner regions is limited by the low 
spatial resolution of the ayailable Hi obseryations. Large disks 
and rings of regularly rotating molecular gas are now known to 



X (arcsec) 

Fig. 5. Brightness distribution, in arbitrary units, along the direction of 
the radio jet of the synchrotron continuum emission (red line) and of 
the CO(2-l) emission (blue line). The dashed vertical lines indicate the 
peaks of the continuum hot spots, showing that the bright CO slightly 
extends beyond the hot spots. 


be relatiyely common structures in early-type galaxies and are 
obseryed in more th an 50% of CO detected early-type galaxies 
(lAlatalo et al.ll201^ . Thus, in this respect, IC 5063 is not an ex¬ 
ception. 

We haye modelled the rota tion of the regu l ar CO distribution 
using the HST photometry of iKulkarni et al.l (Il998h . This pho¬ 
tometry, in the bands FI lOW, F160W and F222M, shows that the 
unobscured light distribution of the inner regions of IC 5063 can 
accurately be described using a de Vaucouleurs’ law. Figure [3 
shows, oyerplotted on a pv diagram, the rotation curye expected 
for the obseryed light distribution where we haye assumed that 
the light traces the mass distribution and that the mass distribu¬ 
tion is spherical. The horizontal scale of the rotation curye is set 
by the obseryed effectiye radius (21 "5 corresponding to 4.3 kpc) 
while we haye set the amplitude of the rotation curye to match 
the rotation yelocities of the yery outer regions in Fig. [3 

4.4. Moiecuiargas masses 

An important property that can be deriyed from the ALMA data, 
are the masses of the yarious kinematical components of the 
molecular gas. The integrated total flux (after primary beam cor¬ 
rection) corresponds to 76.4 + 3.8 Jy km s'. The flux of the gas 
in the region co-spatial with the radio is 38.9 + 1.9 Jy km s 'of 
which 10.0 ± 0.5 Jy kms^'are associated with gas with anoma¬ 
lous yelocities i.e. gas with obseryed yelocities incompatible 
with those of a regularly rotating disk with the rotation curye 
shown in Fig. [3 i.e. gas rotating faster than expected and gas 
that shows apparent counter rotation. The anomalous gas at the 
location of the W lobe has a flux of about 2.7 ± 0.3 Jy km s '. 

In order to conyert these fluxes to masses of molecular gas, 
a numbe r of assumptions are nee ded. As for the APEX data pre¬ 
sented in iMorganti et akl (12013bl) . we estimate the H 2 masses by 
assuming a luminosity ratio L'CO(2-1)/L'CO(1-0) = 1, typical 
of thermalised and optically thick gas, which could be the situa¬ 
tion in the denser and warmer gas associated with an outflow. 

For the kinematically disturbed gas, we haye considered a 
range of yalues for the Xco conyersion factor. This factor is as¬ 
sumed to be between 0.4 x 10^" cm“^ (K km s“')“' (aco = 0.8), 
typical of a thick and dense medium and commonly used for 
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ULIRG, and 0.2 x 10^° (K km (corresponding to 

oco - 0.34, the most conservative case) derived for optically 
thin cases and suggest ed for turbulent gas a ssociate d with winds 
(see th e discussion in iBolatto et alJ (1201 3h and in iGeach et^D 
( 20141) ). Using these two values for the CO-to-H 2 conversion 
factor, we derive masses for the anomalous molecular gas rang¬ 
ing between 1.9 and 4.8 x 10^ M© of which between 0.5 and 
1.3 X 10^ Mq is associated with the gas with the most extreme 
kinematics at the location of the W lobe. The result is, there¬ 
fore, that we confirm that the molecular outflow is more mas¬ 
sive tha^_the_2fi_outflow which was estimated to be 3.6 x 10^ 
M^ dMorganti et alll2007h . We also confirm the large difference 
between the mass of the outflow of cold molecular gas, esti- 
m ated here, and the mass of the warm molecular gas estimated 
bv lTadhunter et al.] (l2014l) . 

From the total integrated flux of the CO(2-l), we derive 
the total mass of the molecular gas of the regularly rotating 
gas outside the region of the radio plasma to be in the range 
Mhj = 1.0 - 1.1 X 10® Mq assuming an Aco conversion fac¬ 
tor equal to 2 x 10^° cm~^ (K km s“')“^ (corresponding to 
oco - 4.6, typical for the Milky Way). We consider the masses 
derived from the ALMA data to be mor e reliable measurement s 
that those obtained from the APEX data (iMorganti et alJ2013bl) . 
Two reasons should be considered. First of all, the blueshifted 
wing of CO was covering most of the available band in the 
APEX data. This made the definition of the continuum and the 
consequen t continuum subtraction very uncertain (as discussed 
already in IMorganti et al.l (l2013bl) l. The large observing band 
of the ALMA data has solved this problem. Furthermore, the 
distribution and kinematics of the anomalous gas as observed 
in the ALMA data has turned out much more complex than 
was assumed when me asuring the gas outflow from APEX by 
IMorganti et alJ(l2013bl) . 

4.5. The radio continuum 

The continuum image at 230 GHz is shown in contours in Fig.|2] 
overlaid to the zoom-in of the central region of the CO(2-l). The 
morphology of the continuum obtained from the ALMA data - 
a triple structure with the core in the middle and two lobes of 
which the western one is the brightest - is very similar to the 
structure observed at lower freque ncies (17 and 24 GHz w ith the 
ATCA) presented and discussed in iMorganti et alJ (l2007l) . Using 
the present and the previous radio continuum observations, we 
have derived the spectral indices of the three components. The 
total fluxes of the three components at 17 and 228 GHz are listed 
in Tahiti] together with the derived spectral indices. Between 
17 GHz and 228 GHz we And a fairly steep spectral index for the 
two lobes (a - 0.87 for the E lobe and a = 0.9 for the W lobe; 
S oc v““) and a flatter spectral index for the core (a = 0.33). Thus 
the values obtained are similar, although slightly steeper, than 
found between the much closer frequencies of the ATCA data. 
As expected, these values support the non-thermal origin of the 
emission with the power law extending to the high frequencies 
of ALMA. The spectral indices also confirm the core and lobe 
classification of the three structures. 

As mentioned above, the coincidence between the radio con¬ 
tinuum and the distribution of the bright molecular gas is clear. 
The CO(2-l) gas wraps around the E lobe and part of the W lobe 
while the brighter CO part appears to overlap (in projection) with 
the W hot-spot. It is interesting that the bright CO region shows a 
strong correspondence with the radio continuum in the part dom¬ 
inated by the most kinematically disturbed gas. As discussed in 


Table 1. Flux densities and spectral indices for the radio three radio 
continuum components in IC 5063. 


Region 

17 GHz 
mJy 

228 GHz 
mJy 

a 

E 

5.6 

0.58 

0.87 

core 

30.8 

12.8 

0.33 

w 

91.2 

6.2 

0.93 



Fig. 6. Position-velocity diagram, as in Fig. [3] with the labels of various 
components identified (described in text). 


Sec. 14.21 this can have influence on both the CO excitation as 
well as production of CR. 

It is also interesting to note that the radial distribution of the 
brightness of CO and radio continuum shows an offset in the 
location of the peaks. The peaks of CO are at slightly larger dis¬ 
tances from the centre compared to the peaks of the radio con¬ 
tinuum on both sides of the radio source, see Fig. 121 We will 
comment on this further below, see Sec. EH 


5. Discussion 

The observed distribution and kinematics of the cold molecular 
gas suggest that the radio plasma jet is driving the bulk of the gas 
outflow. Therefore, we present a simple model which is aimed at 
explaining these characteristics following the scenario prese nted 
in the numerical sim ulations of IWagner & Bicknelll (1201 ll) and 
IWagner et al.l (l2012l) . These simulations describe the effects of a 
newly formed radio jet when moving though a dense clumpy 
medium. A porous medium with dense clumps forces the jet 
to And the path of least resistance, while interacting and grad¬ 
ually dispersing the dense clouds away from the jet axis. In this 
way, clouds can be accelerated to high velocities and over a wide 
range of directions, and along the path of the jet a turbulent co¬ 
coon of expanding gas forms, farther away from the jet axis. 

Before describing the model, we first attempt to And clues 
to properties of the molecular gas in the observed kinematics as 
given in Fig. [3] 
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5. 1. Signatures of jet-ISM interaction in the kinematics of the 
moiecuiar gas in the inner region of IC 5063 

A closer look at Fig. |3] shows that a number of key features can 
be identihed in the molecular gas. These are illustrated in Fig. 
|6] Beyond radii of 2" (i.e. beyond the radio source), we detect 
the large-scale gas disk (labelled D), undisturbed by the radio 
jet. At smaller radii, a bright S-shaped feature is seen. We sug¬ 
gest that this bright feature actually consists of two structures. 
One of these is a linear feature, of about 4 arcsec in size, which 
we labelled R in Fig. |6l When linear features are observed in pu 
diagrams of edge-on galaxies, such as IC 5063, they usually cor¬ 
respond to a rotating, ringlike structure in the galaxy, or to a disk 
with solid-body rotation. Solid-body rotation is not expected in 
the centra l regions of a massive early-type galaxy such as IC 
5063 (e.g. lCappellari et^l2013L see also Fig. |3]l, so it is more 
likely that R corresponds to a ring-like structure. However, Fig. [3 
shows that if R corresponds to a ring, its rotation velocity has to 
be lower than the rotation speed expected from the mass model. 
We suggest that this may indicate that feature R corresponds to a 
rotating ring-like structure with a radius of about 2" but rotating 
at velocities below nominal rotation due to the interaction with 
the outflowing gas. 

At both ends of R, there is bright emission (labelled H) near 
the location of the two radio lobes, but the difference with R is 
that they seem to rotate at the expected rotation speed. Together 
with feature R, these give the S-shape of the bright emission in¬ 
side r = 2". This may indicate that just outside the ring-like 
structure R near the hot spots, there is bright emission, but it is 
rotating at the nominal rotation speed. This is also seen in Fig. |5] 
which shows that the bright molecular gas at the locations of H 
extends slightly beyond the radio hotspots. The excitation con¬ 
ditions of the gas just beyond the hot spots, as seen from the 
centre, may already be changed by radiation from the jet-ISM 
interaction but the kinematics has not been affected yet and the 
gas is still rotating near the nominal rotation velocity. The overall 
kinematics shows that component H only occurs near the radio 
lobes. 

The other feature we identify is the fainter gas inside r - 2" 
with velocities (very) inconsistent with the nominal rotation 
(component O). In particular, part of this component has ’forbid¬ 
den’ velocities which implies that the kinematics of this compo¬ 
nent is not dominated by rotation and that it must involve a fast 
radial component. At the location of the W hot spot the largest 
deviations from rotation are seen (component FO). 

5.2. Modeiing of the kinematics of the gas 

Based on the discussion given above, we attempt to make a geo¬ 
metrical model of the gas kinematics in IC 5063. The basic idea 
of the model is that a plasma jet is expanding in the plane of the 
gas disk of IC 5063 and interacting with this gas disk. This in¬ 
teraction is pushing the gas away from the centre and creating a 
cocoon of laterally expanding gas along the radio jet. 

To model the data cube, we employ a geometrical modelling 
technique often used for modellin g the kinematics of gas disks 
in galaxies (e.g., iJozsa et al]|2007h . i.e. we describe the gas disk 
in IC 5063 by a set of concentric elliptical rings where each ring 
has an orientation, a density, and a rotation- and radial velocity. 
We underline that we attempt to only geometrically describe the 
kinematics of the gas; the model does not include any physics to 
describe, e.g., the excitation of the gas or the actual interaction of 
the jet plasma with the gas disk. For the inclination of the rings, 
we choose 80° and for the position angle 115°. 



Fig. 7. Schematic illustration of the face-on view of the model described 
in Sec. 15.21 with the contours of the continuum emission. 


The face-on layout of the model is shown in Fig.|2]and illus¬ 
trates how IC 5063 would be seen from one of its poles. Along 
the radio jet, the cocoon around the jet drives a lateral expan¬ 
sion (indicated by the black vertical arrows) of the gas of the 
inner gas disk (light brown ellipse in Fig. |7]i away from the jet 
axis. The interaction between the expanding cocoon and the in¬ 
ner gas disk also causes the rotation velocity in this region to be 
well below the nominal rotation, as we seem to observe in Fig. 
|3] This gas would correspond to component O. At the radius 
where this laterally expanding gas meets the undisturbed outer 
disk, a ring-like interface is formed of material with lower rota¬ 
tion velocities, but the expansion is halted. We identify this with 
component R. The gas of the quiescent, large-scale disk outside 
the zone of influence (dark brown Fig. |2]l is rotating according 
to the full rotation expected from the mass distribution (Fig.|3ll. 
The horizontal white arrow Fig.|7]indicates the jet direction. The 
vertical white arrows represent the extreme outflow velocities at 
the location of the W hot spot due to the direct interaction of the 
jet with the ISM (component FO). Just beyond the two hotspots, 
the rotation of the gas in the ISM is following our mass model, 
but has higher brightness due to the nearby hotspot (component 
H). We include this component also in order to reproduce the ob¬ 
servation that the molecular appears to extend somewhat beyond 
the bright radio continuum source (see Fig.|2i. 

We have attempted to reproduce the observed pv diagram of 
IC 5063 by constructing data cubes produced by models of the 
kind described above. Figure [3] shows that the kinematics and 
distribution of the gas is not fully symmetric so any symmetric 
model will have difficulty in fitting the data in detail. Therefore 
we have not attempted to replicate the fine details of the bright¬ 
ness distribution and consequently the quality of a fit is not as¬ 
sessed using a numerical goodness-of-fit criterion such, e.g., a 

Instead, we have focussed on whether a model data cube re¬ 
produces the locations in the pv diagram where gas is detected 
and where it is not. 

We have also attempted to keep the number of free parame¬ 
ters required to make an acceptable match with the data, in the 
sense described in the previous paragraph to a minimum. For ex¬ 
ample, since we access the validity of a model by a qualitative 
comparison with the data, the density/emissivity of the gas of 
components O and R does not have to be a free parameter and 
is assumed to be uniform (i.e. not to depend on position). In this 
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Fig. 8. Position-velocity diagram, as in Fig. [3] with white contours 
showing the model described in Sec. 15:21 


way, the model provides a better test for the assumed kinematics 
and geometry. 

The main two parameters we found to be important to be 
able to describe the data are the velocity of the lateral expansion 
of component O and the reduced rotation velocity of component 
R. Initially, we also varied the orientation of the jet w.r.t. the 
line of sight. We found, however, this parameter to be strongly 
degenerate with the outflow velocity of component O in order to 
reproduce the observed (projected) velocities of component O, 
i.e. for a jet oriented more in the direction of the line of sight, we 
have to increase the lateral expansion velocity in order to obtain 
the same projected velocity while not changing the qualitative 
aspects of the model pv diagrams. We have therefore chosen to 
consider only models for which the jet is moving the plane of the 
sky; such models give a lower limit to the expansion velocity of 
O. For simplicity, we have not included the large-scale disk D in 
our models. 

We have found a good fit to the data using a model as de¬ 
scribed above for the case that the lateral expansion reaches ve¬ 
locities up to 200 km s ' and that the rotation velocity of ring 
R (and of component O) is ~50 km s ' at r = 2" (compared to 
a nominal rotation velocity of 170 kms^' at that radius). Note 
that, since we have assumed the jet to move in the plane of the 
sky, projection effects mean that the 200 km s ' for the lateral 
outflow is a lower limit. We have found it necessary, in order 
to reproduce the large line widths observed for component O, to 
include that this component exhibits a range of expansion veloci¬ 
ties, from 0 km s ' to 200 km s“'. We have chosen to implement 
this by letting the lateral expansion velocity depend linearly on 
the lateral distance from the jet, but other choices, such as, e.g., 
a turbulent medium with a large spread in velocity uncorrelated 
with position, will give very similar results. 

Figure 13 gives the observed pv diagram with superposed on 
it contours based on the above model. Given that we only seek 
qualitative agreement with the data, we only draw those con¬ 
tours that indicate the qualitative properties of the model. Be¬ 
cause we have assumed the spatial brightness distribution to be 
uniform for O and R, the brightness distribution in the pv dia¬ 
gram is the result of kinematics and geometry only, not of any 
density/emissivity structure in the gas disk. It is clear that our 
model reproduces the S-shape of the bright inner feature very 
well. Although the physical brightness of the affected regions O 
and R is uniform in the model, the fact that R forms a ring-like 


structure causes it to project in the pv diagram as a bright ridge. 
The bright regions just outside the hot spots in our model repro¬ 
duce H. Although component O has, in the plane of the galaxy, 
the same brightness as R, it appears with lower brightness in Fig. 
13 because the same amount of emission is spread over a much 
larger velocity range. The match for component O is also very 
acceptable, except for the extreme velocities near the western hot 
spot(F(9). Within the framework of our type of model, this can 
be alleviated by choosing higher lateral expansion velocities at 
these locations. This is justified because at this location the most 
intense jet-ISM interaction is occurring, as evidenced by the very 
bright radio co ntinuum emission from this location. We sug ¬ 
gest, similar to iMomantietaP (ll998|): Oosterloo et al.l d200(lli : 
iMorganti et al.l (l2007h : lTadhunter et al.l ( 20141) . that at the W hot 
spot, the jet is directly interacting with a cloud from the ISM. 

We have considered several variations of the model given 
above. For example, one can assume that the expansion of com¬ 
ponent O is radial from the centre instead of laterally from the 
jet. This would apply for a case where the outflow is driven by 
the central AGN. We found, however, that such models do not 
reproduce the observed pv diagram very well. In particular, in 
such a model, one would expect the observed anomalous veloc¬ 
ities of component O to be largest near the centre and be close 
to zero near the hot spots. This is quite different from observed. 
Another aspect we varied is the shape of the region containing 
component O and of the ring R, i.e. whether it is circular or el¬ 
liptical. Although this also depends on the radial density distri¬ 
bution, we found that it is much easier to explain there is no 
obvious observed radial trend in the brightness of R and O with 
models where R and O have an elliptical shape instead of circu¬ 
lar. This can be roughly understood by considering the extreme 
case of R being a rectangle in which case obviously there is no 
radial dependence introduced by the geometry. 


5.3. Hitting and Pushing: the two effects of the radio jet 

The observed distribution and kinematics of the cold molecular 
gas, as well as the model presented above, support the idea that 
the radio plasma jet is driving the bulk of the gas outflow. 

In this respec t, the model follows the id ea at the base of 
the sim ulations of IWagner & Bicknelll (1201 Ih and IWagner et al.l 
(l2012h . As mentioned above, these simulations describe the ef¬ 
fects of a newly formed radio jet when moving though a dense 
clumpy medium. One should note that, although the clouds do 
initially have densities and temperatures that correspond to a 
molecular phase, the simulations do not include the complete 
cooling of the gas and, therefore, they do not follow exactly the 
same phase of the gas that we observe. However, although a di¬ 
rect comparison of the physical conditions of the gas is not pos¬ 
sible, the kinematics of the gas in the observations and model 
presented here appears to compare remarkably well with the pre¬ 
dictions from the numerical simulations. 

In our case, the clumpy medium is indeed traced by the 
bright CO around the radi o jet. As illustrated by, e.g., Fig. 4 in 
IWagner & Bicknell (1201 ih . the expanding, quasi-spherical en¬ 
ergy bubble is responsible for pushing the clouds outwards. Gas 
outflowing under the effect of the expanding cocoon is expected 
to remain visible all along the jet, in the same way as we have 
observed in IC 5063. 

Interestingly, this causes the jet-ISM interaction to have an 
impact over a much broader region than of the direct cone of the 
jet. This is an important element which makes this kind of inter¬ 
action particularly relevant in the context of feedback effects. 
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The kind of strong impact of a jet we observe would be 
expected when a newly formed radio jet expands into the rich 
ISM surrounding a young radio source. This has indee d been ob¬ 
served in ionized gas and H i in young rad io sources dHolt et al.l 
l2009t[^ih et alJlioOUGereb et al.ll2014bl) . Although we do not 


know whether IC 5063 is indeed a young radio source, the par¬ 
ticular conditions in IC 5063 where the jet is expanding in the 
plane of the large-scale gas disk of the host galaxy, likely makes 
this object an ideal target to test this scenario. 

The velocities of the molecular gas observed with ALMA 
match those obtained by the simulations for jets of compara¬ 
ble power to the one in IC 5063. In addition to this, the ve- 
locity field expected in s uch models, as shown in Fig. 4 of 
IWagner & Bicknelil (1201 Ih . shows that the higher velocities are 
expected at the location where the jet is interacting directly with 
the ISM, as is happening in the W lobe. We claim that at this 
location we are seeing the direct interaction of the jet with a 
large molecular cloud , which is suggested to be present by the 
NICMOS results from Kulkarni et alJ(ll998h . The distribution of 
clouds in the surrounding ISM is likely not symmetrical near the 
two lobes, which could explain why the situation in the E lobe is 
less spectacular. 


5.4. The power source of the outflow 


In the preceding sections we have shown that a plasma jet ex¬ 
panding and interacting in a gas-rich, clumpy medium appears 
to be a good explanation for the morphology and kinematics of 
the gas as observed in our ALMA data. This implies that the ra¬ 
dio plasma jet is the main driver of the gas outflow. However, 
IC 5063 has also a strong (optical) AGN and, therefore, we now 
discuss the phenomena from an energy point of view. 

Using the masses of the outflowing gas given in Sec. 14.41 
we derive the mass outflow rate using M - M/rdyn where M 
is the mass of the outflowing gas and Tdyn a timescale for the 
outflow. Given the complex distribution of the outflowing gas, a 
number of assumptions need to be made. In particular, Tdyn was 
derived using an intermediate distance from the core of 0.3 kpc 
and an average velocity of the outflowing gas of 200 kms^', 
more typical of the gas along the jet than the extreme values at 
the location of the W lobe. This gives Tdyn - 1.5 x 10® yr and we 
obtained M to lie in the range 12 to 30 Mq yr ', depending on 
the assumed CO-to-H 2 conversion factor. 

The kinetic pow er of the outflow ing gas can be estimated 
using the for mula in Holt et all (l2006l) that includes, unlike the 
calculations in ICicone et alJ(l2014l) . both t he radial and turbu lent 
components in the gas motion. Following iHolt et aP (l2006h . we 
have assumed that the relatively large line-width of the outflow¬ 
ing gas reflects a turbulent motion that is present at all locations 
in the outflow region. Thus, the component of turbulent is repre¬ 
sented by the FWHM (full width at half maximum) of the lines: 


£ = 6.34xl0'®y (.L +ergs-'. 

where we have assumed ^out = 200 km s-'and Wturb ~ FWHM ~ 
100 kms-'. These values are giving us a very conservative es¬ 
timate of the kinetic power of the outflow. Combining the Hi 
and CO components, this gives a kinetic power of the outflow of 
-5x10^2 ergs-'. 

We can now compare this with the bolometric luminos¬ 
ity of the AGN and the jet power. The bolometric luminos¬ 
ity of the AGN is estimated to be L®gn “ 2< 10''"' ergs-' 


dMorganti et alJ[2()(y^ . thus the optical AGN is in principle ca¬ 
pable of driving the outflow if the efficiency is higher than 0.6%. 

We estimate the jet pow er u sing the formulae give n in 
IWillot et alJ (flO^ . lwa (1200^ and ICavagnoloetan (1201 Ol) . and 
we obtain values between 5 x 10"'^ erg s-' and of 9 x lO'^ erg s-' 
(using a conservative value of f - 10 for the norm alisation un¬ 
certai nty factor introduced by IWillot et'aP 119991: see also IWul 
I 2 OO 9 I for more details). This indicates that the coupling between 
the energy transmitting to the medium by the jet and the ISM 
should haye high efficiency (between 0.05 and 0.1) in order to 
driye the obseryed outflow. 

It is interesting that both phenomena - AGN and radio jet - 
can power the outflow. Despite of this, the obseryed kinematics 
suggests that the radio jet plays the prominent role in driying 
the outflow. This may indicate that either it has a higher effi¬ 
ciency/coupling, or that perhaps both are at work in the inner 
regions while the radio jet dominates at larger distances. 

It is worth noting t hat IC 5063 was pa rt of the sample in¬ 
cluded in the analysis of ICicone et al.l(l2014l) . Interestingly, how- 
eyer, IC 5063 seems to deyiate from the correlation they found 
between the molecular outflow rate and the AGN luminosity 
(their Fig. 9). This correlation was taken as a signature of the 
optical AGN driying the outflow in the sample of local ULIRGs 
and QSO hosts considered by these authors. The fact that IC 
5063 does not lie on this correlation indicates that the radio jet 
indeed could play a role in IC 5063. 


5.5. Interaction, shocks and cold gas 


Although the model described aboye can explain the kinematics 
of the gas quite well, questions remain about how cold gas can be 
associated with such fast and massiye outflows, since shocks ac¬ 
celerate the gas and may heat it to 10®-10^ K. Molecules sho uld 
dissociate at such temperatures (iHollenbach & McKeelll98^ . In 
the case of shocks prod uced by winds from the accretion disks, 
IZuboyas & Kind (l20I4l) haye proposed that radiatiye cooling is 
fast enough to reform molecules in a large fraction of the out¬ 
flowing material on a short time scale. They assume that, due to 
Rayleigh-Taylor instabilities, the cool gas should collapse into 
clumps, which are then entrained in the hot-diffuse medium. Key 
for the success of such a model is the time scale oyer which the 
clumps can form. 

A similar scenario, albeit with shocks driyen by the radio 
plasma, has been proposed for IC 5063 to expla in the warm 
molecular gas kinematics (iTadhunter et al.l l2014]l and for the 
case o f 4C 12.50 to explain the outflowing Hi ( Morganti et al.l 
l2013ah . In the latter study, a cooling time of the order of lO' 
years was deriyed for the outflowing cloud. Such a time scale is 
comparable to the age of the radio source, suggesting that the 
gas has indeed time to cool during th e process of being expelled. 
Finally, iGuillard et al.l (l2009l l2012h proposed that the turbulent 
gas becomes thermally unstable and cools as it moyes. There¬ 
fore the molecular gas that results from the post-shock cooling 
has the momentum of the tenuous gas that has been accelerated 
to high yelocities. This may explain why the molecular gas is 
characterised by such high yelocities. 

In order to confirm such scenarios, the emissiyity of gas 
strongly disturbed by the passage of a radio jet and associ¬ 
ated shocks needs to be followed with simulations. This has 
been inyestigated in a number of studies, although with par- 
ticul ar emphasis on young stellar and Herbig-Haro objects 

(e.s. IChernin et al.ll994l:lde Gouyeia Dal Pinoll999HRaga et al.l 

l2002h . Inyestigations, using a similar approach, to trace the 
cooler gas by expanding to extragalactic sources the simulations 
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of relativistic jets presented in iRocha da Silval (l2015h . are now 
in progress. Preliminary results (Falceta-Gon 9 alves & de Gou- 
veia Dal Pino priv. comm.) show that a relativistic jet entering a 
denser environment is developing a turbulent cocoon. The con¬ 
sequent fast cooling, associated with the turbulence (plus en¬ 
trainment), results in a mix of gas phases (from hot to cold) as 
is observed. Most importantly, the simulations show, to first or¬ 
der, that the cooler gas is also characterized by high velocities. 
Thus, despite the simplicity of the assumptions made (i.e. ho¬ 
mogeneous ambient medium), the shocked ambient material de¬ 
posited into the cocoon around the jet includes a cool component 
(in addition to the warm and hot) with large turbulent velocities, 
in agreement with the observations. 

The results of our observations can also be used to investi¬ 
gate alternative scenarios. For example, slow entrainment by the 
jet that would not dissociate the molecules. However, if the warm 
H 2 is also heated during this entrainment, the mass discrepancy 
between the warm H 2 (Mhj ~ 8.2 x 10 ^ Mq, iTadhunter et alJ 
(120141) ') and cooler component derived from the CO(2-l) obser¬ 
vations presented here, is difficult to explain. Alternatively, if the 
heating of the clouds (and of the warm molecular gas) is done 
by the AGN radiation field (i.e. in a similar way as seen in ra- 
diatively excited molecular clouds, where the warm H 2 forms a 
skin around the colder molecular material deeper in the clouds). 
However, our new CO(2-l) observations make this possibility 
unlikely, because it is difficult to explain how more than 10 ^ 
Mq of molecular gas (equivalent to an extremely massive Giant 
Molecular Cloud) can be accelerated to the observed high veloc¬ 
ities without heating the molecules and destroying the clouds. 

Thus, we conclude that the most likely explanation for the 
outflowing molecular gas in IC 5063 is rapid post-shock cooling. 


6. Conclusions 

The ALMA data of IC 5063 presented here have provided a 
unique, high-resolution 3-D view of the molecular gas in one 
of the best examples of an on-going interaction between the ra¬ 
dio jet and the ISM. Because of its proximity, and its geometry, 
this object is ideal to study the impact of radio jets on the sur¬ 
rounding medium. A wealth of data, including the warm molec¬ 
ular gas, HI and ionized gas, are already available for this object. 
However, the ALMA data have provided a much more complete 
and complex picture and have allowed to trace this interaction 
in much more detail (down to scales of 0.5 arcsec, i.e. 100 pc) 
compared to what is available for the other phases of the gas. 
The results presented here further confirm the multi-phase na¬ 
ture of AGN-driven outflows and the significant mass of the cold 
gaseous component associated to it. 

From these data, we have found kinematically disturbed gas 
distributed all along the radio emission, with velocities strongly 
deviating from the regularly rotating gas. This illustrates that 
gas outflows are produced not only where the jet directly in¬ 
teracts with the ISM, but that the effects occur over a much 
larger volume therefore making its impact larger, as afreadv 
shown by numerica l simulations (e.g.. l'Wagner & Bicknellll2()Tll : 
I'Wagner et al.ll20l3) . 

The region showing the disturbed kinematics of the gas is co- 
spatial with the radio plasma and it also appears much brighter 
than the gas of the larger-scale, regularly rotating component of 
molecular gas. This inner, brighter region of CO(2-l) is elon¬ 
gated along the radio axis and wrapping around the radio contin¬ 
uum emission. The brightest CO(2-l) emission is co-spatial with 
the bright, western radio lobe and has the largest velocities de¬ 
viating from regular rotation. A simple model assuming lateral 


expansion of the gas caused by the cocoon around the radio jet, 
describes the observed kinematics remarkably well. 

Taking all this together, we argue that the distribution and 
kinematics of the molecular gas support the idea of the jet inter¬ 
acting with a clumpy medium and pushing it outwards. In addi¬ 
tion to the outflowing gas, this results in the creation of a ring of 
gas, which represents the interface between the region affected 
by the radio plasma and the undisturbed disk at larger radius. At 
least at the location of the western hot spot, the jet is directly in¬ 
teracting with a gas cloud, causing the very large outflow veloc¬ 
ities observed. Suc h a picture is fully consiste nt w ith the numer¬ 
ical sim ulations of lWagner & Bicknelll (1201 ll) and I'Wagner et al.l 
(I 2 OI 2 I) of which all main components are clearly identifiable in 
the data. 


Even assuming the most conservative value of the conversion 
CO-to-H 2 , we find that the mass of the outflowing gas is between 
1.8 and 4.5 x 10^ Mq, thus much larger than the outflow observed 
in any other gas phase. The mass outflow rates we derive are in 
the range 12 - 30 Mq yr“^ 

Interestingly, looking at the energetics, it appears that both 
the bolometric luminosity of the AGN, as well as the jet power, 
are large enough to drive the outflow. The finding, based on the 
kinematics, that the radio plasma is playing the dominant role 
(at least at radii larger than 100 pc), could suggest a higher effi¬ 
ciency for this mechanism. 

Particularly relevant from this study is the conclusion that 
the effect of the radio plasma can be significant, also in objects, 
such as IC 5063, that are often considered radio quiet. A number 
of cases where this situ ation may occur have been recently iden - 
tified (e.g. NGC 1266 dAlatalo et akllMlt iNvland et aljf2013l) : 
NGC 1433 dCombes et alJl2013[): .11 430-H1339 - nicknamed the 
“Teacup AG” - dHarrison et al.l 20141) '). However, direct evidence 
that the radio plasma is at the origin of these outflows of cold 
gas is, so far, not available in all these cases. Another case is 
M51 where signatures of the interaction between a weak radio 
jet and the ISM have been found in the molecular gas. A high- 
density region traced by HCN has been detected along the weak 
radio jet of this nearby object. This region shows a gradient in 
the HCN/C O ratio, with decreasing values for larger distances 
from the jet (iMatsushita et aljr20L5l) . This points to the effect of 
the jet pushing and compressing the surrounding medium. If the 
possibility of weak radio jets driving massive gas outflows, due 
to their high coupling efficiency, is confirmed for more objects, 
it would make even “wimpy” radio jets relevant for the fate of 
gas in galaxies. Because such weak radio jets are more common 
than the powerful ones, this can have implications for the role of 
feedback in galaxy evolution. 


Another interesting point in the context of feedback effects is 
the fate of the kinematically disturbed gas. The observed outflow 
velocities of most of the disturbed gas are only ~200 km s '. Un¬ 
der such conditions, the gas will not leave the galaxy, but most 
likely will just be relocated, in a process comparable to a “galac¬ 
tic fountain”. Of the gas with disturbed kinematics, only at most 
a very small fraction may have velocities high enough to leave 
the galaxy, a simila r situation as observed in, e.g., NGC 1266 
(lAlatalo et al.l2015l) . However, the injection of the large amounts 
of energy in the gas disk may have a similar effect as increasing 
the turbulence of the gas, a process that is also considered to in - 
hibit star formation (see the case of 3C326. lGuillard et al.l2015l) . 


In summary, the ALMA data presented here for IC 5063, 
show how spatially resolved observations are key to being able 
to determine which mechanism, radiation or jet, drives the ob¬ 
served kinetic feedback. Hopefully, expanding this kind of work 


Article number, page 10 of 1 11 1 





































Morganti et al.: The fast molecular outflow in the Seyfert galaxy IC 5063 as seen by ALMA 


to more objects will finally allow us to quantify the impact of an 
AGN on the ISM of the host galaxy. 
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